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High field magnetization, field- dependent specific heat measurements, and zero field inelastic mag- 
netic neutron scattering have been used to explore the magnetic properties of copper pyrazine 
dinitrate (Cu(C4li4N2)(N03)2). The material is an ideal one-dimensional spin-1/2 Heisenberg anti- 
ferromagnet with nearest neighbor exchange constant J = 0.90(1) meV and chains extending along 
the orthorhombic a-direction. As opposed to previously studied molecular-based spin-1/2 magnetic 
systems, copper pyrazine dinitrate remains gapless and paramagnetic for g^sH/ J at least up to 1.4 
and for ksT / J at least down to 0.03. This makes the material an excellent model system for ex- 
ploring the T — critical line which is expected in the H — T phase diagram of the one-dimensional 
spin-1/2 Heisenberg antiferromagnet. As a first example of such a study we present accurate mea- 
surements of the Sommerfeld constant of the spinon gas versus g^sH/J < 1.4 which reveal a 
decrease of the average spinon velocity by 32% in that field range. The results are in excellent 
agreement with numerical calculations based on the Bethe ansatz with no adjustable parameters. 
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I. INTRODUCTION 

Quantum many body systems which support qualita- 
tively different ground states as a function of Hamilto- 
nian parameters have the potential to display so-called 
quantum critical phenomena when these parameters are 
tuned in the vicinity of zero temperature phase transi- 
tions. Many unusual physical phenomena such as 
heavy fermion behavior, || high temperature supercon- 
ductivity, [^,^| and non-fermi-liquid behavior || are now 
thought to be manifestations of proximity to such quan- 
tum critical points. To understand the intrinsic proper- 
ties of quantum criticality it is important to explore its 
phenomenology in simple and well controlled model sys- 
tems. One of the simplest quantum critical many body 
systems is a linear chain of antiferromagnetically coupled 
spin-1/2 objects. Quantum criticality in this system is 
particularly interesting because it is possible to continu- 
ously tune the critical exponents by the application of a 
magnetic field. 

Unlike both conventional ordered magnets and gapped 
quantum spin chains which generally have sharp, low- 
energy dispersive modes, the excitation spectrum of 
the S=l/2 antiferromagnetic (AFM) chain is dominated 
by the so-called "spinon continuum," [Hjjn]] which is 
sketched in Fig. [I]. Some properties of the S=l/2 AFM 
chain may be computed exactly using the Bethe ansatz, 
jl2| but it is also very illuminating to map the spin chain 
onto a ID system of interacting fermions. [ ft3| This gives 
important insight into the spinon continuum, which may 



be viewed as the particle-hole continuum of the fermion 
model. This underlying fermion character of the S=l/2 
chain also has consequences for thermodynamic proper- 
ties such as the magnetization and specific heat. There 
are a growing number of magnetic materials in which the 
predictions of these models have been tested experimen- 
tally. 
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FIG. 1. Schematic view of the bounds of the spinon con- 
tinuum of the S=l/2 AFM chain, drawn for J = 0.90 meV as 
appropriate for CuPzN. Solid dots show fixed values of q for 
the data in Fig. ^. Horizontal and vertical lines show the scan 
directions for the data in Figs. ^ and ^. The ellipse indicates 
the FWHM instrumental resolution for Figs. ^ and ^. 

In the presence of a magnetic field H , an isotropic spin 
chain is described by the Hamiltonian 



(1) 



The field causes a Zeeman splitting of the half-filled, 
doubly-degenerate spinon band, which introduces addi- 
tional low-energy spanning vectors, and the spin chain 
has been predicted (ll]] to develop new soft modes at 
wave- vectors which are incommensurate with the lattice, 
but which connect the field-dependent Fermi points. The 
dynamic critical exponents which control the tempera- 
ture dependence of spin fluctuations close to these soft 
points are predicted to vary continuously with applied 
field. In a recent experiment on copper benzoate we have 
obtained the first direct experimental evidence for the 
existence of these low-energy modes. |Q However, while 
this work showed directly the field-induced incommensu- 
rability of the S=l/2 chain, it also showed a field- induced 
gap in copper benzoate. Oshikawa and Affleck showed 
p9| that this gap arises from the presence of two Cu 
sites per unit cell along the copper benzoate chains. The 
combination of a staggered g-tensor and Dzyaloshinskii- 
Moriya interactions induces an effective staggered field 
when a uniform field is applied and this in turn creates 
a gap in the excitation spectrum and other interesting 
changes in the dynamic spin correlations. fl9| , |20|| 

The result of this is that copper benzoate does not 
have a quantum critical line in the H — T plane be- 
cause a field immediately drives the system away from 
criticality. In this paper we describe measurements on 
a long neglected one-dimensional spin-1/2 AFM, cop- 
per pyrazine dinitrate (CuPzN), which remains gapless 

for g^isH ~ J and fc^T <C J, and thus is very well 
suited for experimental studies of quantum criticality. 
CuPzN is a S=l/2 chain with intra-chain exchange con- 
stant J « 0.91 meV (J/fcs ~ 10.6 K), as determined 
from previous zero-field susceptibility and specific heat 
studies. |^l],^2| It is highly one-dimensional; 3D order- 
ing having not been observed in zero field for T > 0.1 
K, which implies [^3| that the ratio of inter-chain 
to intra-chain coupling is J' /J < 10~ 4 . The magnitude 
of the exchange constant J in CuPzN makes the mate- 
rial an attractive candidate for both thermodynamic and 
neutron scattering studies in magnetic fields, as large ef- 
fective fields gfisH/J can be obtained. We report high 
field magnetization, field-dependent specific heat, and 
zero-field inelastic neutron scattering measurements that 
show that CuPzN is a nearly ideal example of an isotropic 
Heisenberg AFM chain. In addition, our high field spe- 
cific heat data provide direct experimental evidence for 
the suppression with field of the spinon velocity in the 
one-dimensional spin-1/2 AFM. 



II. COPPER PYRAZINE DINITRATE 

CuPzN, Cu(C 4 H 4 N2)(N03)2, is orthorhombic, with 
space group Pmna and room temperature lattice con- 
stants a = 6.712 A, b = 5.142 A, and c = 11.732 A. pi 



The copper ions form chains along the [100] direction, 
with one copper per unit cell along the chain. The Cu 
ions along the chain are coupled magnetically through 
the pyrazine molecules, as shown in Fig. |^. The room 
temperature g-factors are g a — 2.05, gb — 2.27, and 
g c = 2.07. || 




FIG. 2. Crystal structure of copper pyrazine dinitrate, 
Cu(C4H 4 N 2 )(N03)2, showing how Cu 2+ ions (hatched 
spheres) are linked through pyrazine rings to form ID chains. 
The chain axis (a) is vertical on the page, with the b axis 
nearly horizontal. The long bonds between the 02 atoms 
and the copper sites are represented by thin lines (intrachain 
contacts) and dashed lines (interchain contacts). Atomic po- 
sitions are given in Table 1, and bond lengths and angles are 
given in Table 2. 

We have recently re-determined the cry stal structure 
of CuPzN with X-rays at T = 158 K. |§ There are no 
major differences from the room temperature structure 
as originally determined. [p4[ The cell parameters show 
very slight contraction upon cooling along the chain axis 
a (6.712 A to 6.697 A ), and somewhat more change 
along the two axes perpendicular to the chains (b axis: 
5.142 A to 5.112 A; c axis: 11.732 A to 11.624 A). With 
neutrons, we have determined the low-temperature (T < 
4 K) lattice constants to be a — 6.6870 A, and b = 5.0983 
A. Throughout this paper we index wave vector transfer 



in the corresponding reciprocal lattice: Q = ha* + kh* + 
lc*, where a* = (2n/a)k etc. We also define wave vector 
transfer along the chain as q = Q • a = 2nh. 

The atomic coordinates and some of the important 
bond lengths and bond angles measured at T — 158 K are 
listed in Tables | and O. Each copper ion is coordinated 
to two nitrogen atoms from the adjacent pyrazine rings 
(Cu-Nl = 1.973 A) and two oxygen atoms (01) from 
two different nitrate ions. The Cu-01 distance equals 
2.004 A and is represented by a thick solid line in Fig. [|. 
The copper atoms and nitrate groups lie in mirror planes 
perpendicular to the two-fold a axis. There is a semi- 
coordinate bond between each copper atom and a second 
oxygen (02) from each nitrate group, represented by a 
thin line in Fig. f| The Cu-02 distance is 2.478 A, which 
is much greater than the 1.9 A for copper oxides with oc- 
tahedral coordination, and the angle between the Cu-02 
bond and the normal to the Cul-Nl-01 plane is 30.4°, 
which is also unfavorable for superexchange interactions. 
The 02 atoms also bridge the copper atoms to the cop- 
per sites on adjacent chains by a long pathway (3.264 A, 
dashed lines). The axis connecting the Cu site to the 
02' site on the adjacent chain is tilted only 21.6° from 
the normal to the local coordination plane. This long 
interaction may help stabilize the formation of the three 
dimensional lattice, but there is no evidence that it pro- 
vides a superexchange pathway. The basal plane of the 
copper site is determined by the four short bonds. The 
strength of the superexchange interactions through the 
pyrazine rings will be a function of the angle between the 
normal to the basal plane and the normal to the (planar) 
pyrazine molecules which is 51°. [E7| 
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TABLE I. Fractional atomic coordinates and equivalent 
isotropic thermal parameters for CuPzN at T = 158 K. 



Bond Length 


Bond Angle 


Cul-Nl = 1.974(3) A 
Cul-Ol = 2.004(3) A 
Cul-02 = 2.478(3) A 


Nl-Cul-Ol = 90° 
Ol-Cul-02 = 56.88(9)° 
Nl-Cul-02 = 90° 



TABLE II. Selected bond distances and angles for CuPzN 
at T = 158 K 



III. EXPERIMENTAL TECHNIQUES 

Crystals of CuPzN were grown by slow evaporation of 
aqueous solutions of Cu(ll) nitrate and the corresponding 
pyrazine (either protonated or deuterated). Large single 
crystals formed after several months. Faster evaporation 
results in the formation of many small, needle-shaped 
crystals. 

Specific heat measurements on small, protonated crys- 
tals with typical mass m — 5 mg were made in fields 
up to fi H = 9 Tesla for 0.1 K < T < 10 K using the 
relaxation method. |^8|,^] Magnetization measurements 
on both protonated and deuterated powder samples were 
carried out at the National High Field Magnet Labo- 
ratory in magnetic fields up to hqH — 30 T using an 
EG&G vibrating sample magnetometer which had been 
calibrated against a high-purity nickel sample. 

For inelastic neutron scattering, a sample was used 
which consisted of two deuterated crystals with total 
mass m = 72.8 mg, co-aligned to within 0.4° in the (hkO) 
scattering plane. The measurements were performed on 
the SPINS cold neutron triple axis spectrometer at NIST. 
For measurements requiring a nearly isotropic wave vec- 
tor resolution function we used a conventional parallel 
beam configuration. In this configuration the projection 
of the instrumental resolution ellipsoid on the scattering 
plane is an ellipse with principal directions approximately 
parallel and perpendicular to the wave vector transfer Q, 
and a Full Width at Half Maximum (FWHM) along these 
directions of typically AQy = 0.02 A" 1 and AQ±_ = 0.054 
A -1 respectively. For measurements requiring good Q 
resolution only along a single direction we used a hori- 
zontally focusing analyzer. With this configuration the 
principal axes of the projection of the resolution ellipsoid 
on the scattering plane are approximately parallel and 
perpendicular to the direction of the scattered beam k f , 
and the FWHM were AQy = 0.045 A -1 , and AQ ± = 
0.32 A -1 respectively. By orienting the spin chain axis 
parallel to k/ the instrument effectively integrates over 
wave vector transfer perpendicular to the chains. This 
led to a gain in intensity by a factor of 5.3 over the con- 
ventional parallel beam configuration without significant 
loss in wave vector resolution along the chain axis. The 
projection of the resolution ellipsoid on the energy trans- 
fer axis had FWHM ATiuj = 0.15 meV in both configu- 
rations. The half-value contour of the projection of the 
resolution function of the focusing analyzer configuration 
onto the Tiuj-q plane is shown in Fig. [jj We used 20 cm 
of beryllium oxide at T = 77 K after the sample to re- 
ject neutrons with energies above the fixed final energy 
Ef = 3.7 meV from the detection system, and 10 cm of 
beryllium at 50 K before the sample to suppress higher 
order contamination of the incident beam for incident 
energies £7, = Et + fiio < 5.1 meV. The detector count 
rate was measured in units of the count rate of a monitor 
with sensitivity proportional to the incident wave length. 
With this standard technique the raw count rates ex- 



shifts to lower T. The solid lines are the results of a fit to 
a model based on exact diagonalization of short chains, 
as described below. The dashed line is the estimated 
phonon contribution determined from the fit. 



D 



□ 




T(K) 

FIG. 4. Specific heat of CuPzN vs T at constant magnetic 
field fi H = 0, 3 T, 5 T, 7 T and 9 T || b. Solid lines are 
a fit to an exact diagonalization model as described in text. 
Dashed line is the phonon contribution determined from the 
fit. 
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Figure |^ shows an expanded view of the data in Fig. |] 
for T < 1 K, plotted as C/T vs. T 2 . In addition, the 
specific heat at ^qH = 9 T || a is shown. A predom- 
inantly linear in T behavior is observed. The positive 
slope visible at /jqH — 9 T is a consequence of the field- 
induced shift of spectral weight to lower temperature. An 
increase in the linear part of C (T) with increasing field is 
also seen. This indicates a re-normalization of the spinon 
velocity with field, as will be discussed below. 



FIG. 3. Relative magnetization M/M sa t of CuPzN at 
T = 1.8 K and 4.2 K. Dot-dashed line: calculated T = 
magnetization (Ref. 0) of an S=l/2 AFM chain with 
J/k B = 10.3 K. Dashed line: fit of T = 1.82 K data to fi- 
nite chain model described in text. The fit of the T = 4.2 K 
curve is indistinguishable from the data. Inset: fit residuals 
SM — M — M fit . Solid line: T = 1.8 K; dashed line: T = 4.2 
K. 



B. Specific Heat 

Figure ^ shows the temperature dependence of the to- 
tal specific heat of CuPzN for a number of magnetic fields 
H || b. As H is increased, the broad feature observed 
in zero-field is suppressed, and the maximum gradually 




FIG. 5. Specific heat of CuPzN below T = 1 K, plotted as 
C/T vs T 2 for at constant magnetic field fi Q H = 0, 3 T, 5 T, 
7 T and 9 T || 6, and fj, H = 9 T || a. 



C. Neutron Scattering 

While bulk measurements have established that 
CuPzN is a quasi-onc-dimcnsional magnet they cannot 
distinguish along which directions spins are coupled and 
along which directions they are not. For CuPzN both 
the a and b directions are candidate chain directions. To 
distinguish between these possibilities we measured the 
wave vector dependence of the inelastic magnetic neu- 
tron scattering cross section in the (hkO) plane at low 
temperatures and for low values of energy transfer hui. 
Fig. ||(a) shows three scans along the b* direction col- 
lected at T — 0.3 K and hui = 0.29 meV, and for three 
values of wave vector transfer along the perpendicular 
a* axis (indicated as dots in Fig. [l]). None of the scans 
shows any statistically significant modulation. There is, 
however, clearly more intensity for Q ■ a* = tt than for 
the other nearby values of this component of wave vector 
transfer. This is brought out by Fig. |6|(b) which shows 
the difference between the Q • a* = tt and the Q • a* ^ tt 
data. No modulation beyond that associated with the 
single-ion squared magnetic form factor (solid line) is ob- 
served, which indicates the absence of spin correlations 
along the b* direction on a time scale 1/ui — 2 ps for 
ksT <C hu. This observation allows us to use the hori- 
zontally focusing analyzer to gain sensitivity by probing 
the b*-integrated intensity as we explore the q = Q • a 
dependence of <S(Q,w). 
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FIG. 6. Neutron scattering intensity vs wave vector trans- 
fer fcb* perpendicular to the chains at hui — 0.29 meV at 
T = 0.3 K. (a) raw data at h = 0.5, 0.38, and 0.62 (q = tt, 
0.767T, and 1.24.7T.) (b) Difference AJ between h = 0.5 scan 
and the average of the h = 0.38 and h = 0.62 scans. The lack 
of variation in A7 with k shows the one-dimensionality of the 
magnetic scattering in CuPzN. 

Data collected with the focusing analyzer are shown 
in Figs. and ||. Figure ^ shows the inelastic scattering 
intensity versus wave-vector transfer q along the chains 
at T = 0.3 K for three fixed values of energy transfer 
hu = 0.29 meV, 0.54 meV, and 0.79 meV. Note that in 



all these scans wave vector transfer along the b* direction 
was varied so as to maintain kf || a* as required to main- 
tain good wave vector resolution along the spin chain in 
the focusing analyzer configuration. At hui = 0.29 meV 
a single peak is observed at q = tt, consistent with the 
data in Fig. [| When the temperature is raised to T = 20 
K this peak disappears, as shown by the solid symbols 
in Fig. 0(c). This indicates that the inelastic scattering 
observed is indeed magnetic in origin. As hu is increased 
at T — 0.3 K two peaks become visible, demonstrating 
the presence of dispersion along the chain direction. At 
hui = 0.79 meV, the two peaks are clearly resolved, but it 
is also apparent that the scattering cross section between 
the peaks is finite which indicates a continuum contribu- 
tion. The locations of the maxima in these constant-?Lcj 
scans coincide with the lower bound of the spinon con- 
tinuum anticipated for CuPzN on the basis of bulk mea- 
surements (see Fig. [l]). Figure |^ shows the scattering 
intensity vs. hu at q — tt/2 (ie. along the vertical line in 
Fig. |). The abrupt increase in intensity at hu ~ 1.3 meV 
and the subsequent gradual decay at higher energies are 
again characteristic of the bounded spinon continuum. 
The solid lines in Figs. ^ and || were obtained from a fit 
to a model of the two-spinon contributions to S(q,u>) to 
be described below. 
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FIG. 7. Neutron scattering intensity vs wave vector trans- 
fer along the chains. Open symbols: T = 0.3 K. Filled sym- 
bols: T = 20 K. Right-hand axis gives normalized intensity 
scale as described in text. Solid lines are fit to a T = model 
farS(«,w) (Ref. @). 
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FIG. 8. Neutron scattering intensity vs energy transfer at 
q = 7r/2 and T = 0.3 K. Right-hand axis gives normalized 
intensity scale as described in text. The solid line is fit to a 
T = model for S(q,cu) (Ref. @). 



V. DISCUSSION 
A. Magnetization 

As is well known, many of the finite temperature ther- 
modynamic properties of quantum spin chains can be 
calculated quite accurately from exact diagonalization of 
short chains. (3^] Comparison of data to calculations of 
this type is a very effective way of ascertaining how well 
a material is described by a particular model, and also of 
determining the exchange constants in that model. We 
now compare our results for the magnetization and spe- 
cific heat of CuPzN to exact diagonalizations of linear 
chains containing up to N = 16 spins-1/2. 

For the magnetization, we find surprisingly good agree- 
ment between our data and the calculated result for 
N = 16. We carried out a simultaneous fit to both 
data sets shown in Fig. | (T = 1.82 and 4.2 K), us- 
ing the measured average g-factor g avg = 2.11. This fit 
gave J/ks — 10.3(1) K. |33[ The results of the fit are 
superimposed on the data in Fig. ^| but can only be dis- 
tinguished from the data near [IqH = 15 T at T = 1.82 
K. The residuals SM = M — Mfu are shown in the in- 
set of Fig. g[ These show the deviation of model from 
data to be significantly less than 1% over most of the 
field range. As a measure of the error in this fit associ- 
ated with the finite system size used in the calculation, 
we have studied the convergence of M(H, T, N) with in- 
creasing N. At T — 1.82 K (k B T/J = 0.17), the dif- 
ference [M (N = 16) - M(N = 15)]/M sot is less than 
0.003 for all H. At T = 4.2 K, this difference is less than 
2 x 10~ 5 . We conclude that our exact diagonalization 
results are representative of infinite length chains in the 
parameter range studied and that the magnetization of 
CuPzN down to ksT/J — 0.17 is perfectly described by 
the linear chain model. 



B. Specific Heat 

At the high end of the temperature range accessed in 
our measurements, the specific heat may also be calcu- 
lated accurately from exact diagonalization. Following 
Bonner and Fisher, p3] we compute the specific heat at 
constant field Ch(J/T) as the average of that obtained 
from chains with 15 and 16 spins. Finite-size scaling 
suggests that this procedure gives quite accurate results 
down to temperatures ksT/J ss 0.2. All the data for 
T > 2K in Fig. | were fit simultaneously to 



C(T) = AC H (J,T) + BT 3 



(2) 



where the second term accounts for the lattice contri- 
bution. As shown by the solid lines in Fig. ||, this 
model accounts for the data very well, with A = 1.05(1), 
J/k B = 10.6(1) K, and B = 2.4(1) mJ/molc-K 4 . We 
ascribe the deviation of A from 1 to a systematic error 
in the normalization of our specific heat data. This is 
corroborated by a computation of the magnetic entropy 
from the specific heat data, which comes out to be 4% 
higher than the expected value of R In 2. [p9| 

The linear dispersion relation of the fermionic spinons 
at low energy implies that the low-temperature specific 
heat of the S=l/2 AFM chain should be linear in T. It 
is given by [ ^4|j35| 



C H {T)= 1H T = \R^- V 
3 v s {H) 



(3) 



where v s (H) is the field-dependent spinon velocity. In 
zero field, v s — 7rJ/2, and the specific heat becomes C = 
(2/3)RksT/ J, (3^^] a result which has been confirmed 
experimentally in both copper benzoate |Q and CuPzN. 
fl22fl The field-dependence of v s (H) can be determined 
from the Bethe ansatz |7],|38| and used to predict Ch {T) . 
It is calculated as 



E 



27rcr(A) ' 



(4) 



where E and a are determined by the following set of 
integral equations 



1 - ^ / K(r)- rfW)H, 

27r J -A 



7T(1 +7/ 2 ) 



27r J -A 
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dK{r] - A) 



drj 



-A 



E = 



4A 



(1 + A 2 ) 2 



e (r))p(r])dr], 



with K(rj) — 4/(4 + rf ). The boundary of integration 
A is related to the magnetic field H by g^sH/J = 
27rcr(A)/£(A). In practice, we first fixed the value of A, 



and then solved the integral equations by numerical in- 
tegration and iteration to obtain v s and H. Accuracy of 
more than six digits can be achieved within minutes on a 
standard workstation. We performed the calculation for 
many values of A to give v s as a smooth function of i?, 
which is shown in the inset to Fig. ^. This determines 
the theoretical curve of the coefficient 7(if ) with no ad- 
justable parameters, once J and the g-value are fixed. 

Figure [| shows the measured field dependence of 7^ = 
Ch/T determined from the data in Fig. ||, together with 
the calculated result. These are plotted against reduced 
field H* = gfi B H/J with J/k B = 10.6 K, and g a = 2.05 
and gi, — 2.27. [g5| Note that for this comparison a 
correction factor 1/A = 0.95 was applied to the spe- 
cific heat data to account for the normalization error 
discovered when comparing to exact diagonalization cal- 
culations. The agreement between model and data is 
excellent and the data provides direct experimental evi- 
dence for the field-dependence of the spinon velocity in 
the S=l/2 AFM chain. 
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FIG. 9. Field dependence of the Sommerfeld constant 
(7 = liniT^o C/T)derived from low-T specific heat data, com- 
pared to Bethe ansatz calculation (solid line). The data was 
scaled by a factor 0.95 to account for a systematic error in 
normalization which was derived from the comparison to ex- 
act diagonalization data shown in Fig. ^. Inset: calculated 
field-dependence of the spinon velocity v B (H). 



C. Neutron Scattering 

The magnetic contribution to the measured inelastic 
neutron scattering signal is proportional to the normal- 
ized magnetic scattering intensity 

7(Q,c) = ||F(Q)| 2 ^(5 a/3 -QaQ^)5 a/3 (Q,w). (5) 

a/3 

In this expression F(Q) is the magnetic form factor of 
Cu 2+ , fH and S af3 (Q,uj) is the convolution of the dy- 
namic spin correlation function, p0[ 
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2nh 



die^lg<5|(t)«S|,(0)>e- 

RR' 



■(R 



(6) 



with the normalized instrumental resolution function. 
pl[ We determined J(Q,w) in absolute units by com- 
parison of the magnetic scattering intensity to the inco- 
herent scattering intensity Iy from a vanadium reference 
sample: 



1 



Hy 



J dfnoly n Cu r 



I mo g(Q,w). (7) 



Here ny and ncu are the number of moles of vanadium 
and copper in the reference and sample, respectively, 
hi — 6.36 fm is the incoherent scattering length for vana- 
dium, p0| ro = 5.38 fm, and I m ag(Q, is the portion of 
the measured signal attributable to magnetic scattering. 
The conversion of the measured count rates to I(0,oj) 
is shown on the right-hand axis in Figs. and |8[ As 
a complete, independent measure of the non-magnetic 
background was not made, the zero of the I(Q,iv) scale 
was set from the backgrounds determined in fitting the 
data. 

Miiller et al. |Tl| l proposed the following approximate 
expression for the T — dynamic spin correlation func- 
tion for the spin- 1/2 Heisenberg AFM chain, 



1 



A 



2tt 



V(M 2 -*f(g) 



= Q(hto - e 1 {q))e{e 2 (q) - hcu) 



(8) 



In this expression Q(x) is the step function, while £i(q) = 
(ttJ/2)\ sinq| and £2(9) = itJ sin q/2 are the lower and up- 
per bounds of the spinon continuum. Recent publications 
present exact calculations of the contribution to S(q, lo) 
from n-spinon excitations. Jl(],|^,|]| The numerical dif- 
ferences between the exact expressions for S(q, uj) and 
the simpler approximate form are small, and are unim- 
portant for the analysis of our relatively low statistics 
scattering data for CuPzN. For the purpose of determin- 
ing whether the scattering data can be described by the 
spin chain model which accounted for our bulk data we 
have therefore used Eq. ||. This expression was convolved 
with the experimental resolution function and the overall 
scale factor, A, the exchange constant, J, and a constant 
background for each scan were varied to achieve the best 
possible fit to the scattering data. The results of this 
fit with A = 3.5(1) and J = 0.90(1) meV are shown as 
solid lines in Figs. |?] and |[ and may be seen to provide 
an excellent description of the data. The value of J de- 
rived from neutron scattering is indistinguishable from 
that derived from specific heat data. 



VI. CONCLUSION 

The combination of magnetization, specific heat and 
inelastic neutron scattering reported here show that 
CuPzN is extremely well described by the simple 
isotropic Heisenberg model of the S=l/2 AFM spin chain 
with J = 0.90(1) mcV. The neutron scattering mea- 
surements have confirmed that the crystalline a-axis is 
the ID chain axis, as had been predicted from struc- 
tural analysis. The field-dependence of the specific heat 
shows that CuPzN remains gapless in a field and that the 
spinon velocity is reduced in agreement with predictions 
based on calculations using the Bethe ansatz technique. 
Taken as a whole, these results establish that CuPzN 
is an excellent model system for the one-dimensional 
spin-1/2 Heisenberg antiferromagnet even in large mag- 
netic fields (gjiBH/ J ~ 1.4) and at low temperatures 

(ksT/J ~ 0.03). This material is therefore an excellent 
candidate for further studies of the unique field depen- 
dent quantum critical properties of the spin-1/2 Heisen- 
berg antiferromagnet. 
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